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Abstract 

A study of the pd — > pdr] reaction in the energy range where the recent data from Uppsala 
are available, is done in the two-step model of r] production including the final state interaction. 
The rj — d final state interaction is incorporated through the solution of the Lippmann Schwinger 
equation using an elastic scattering matrix element, T^^^^^, which is required to be half off-shell. 
It is written in a factorized form, with an off-shell form factor multiplying an on-shell part given 
by an effective range expansion up to the fourth power in momentum. The parameters of this 
expansion have been taken from an existing recent relativistic Faddeev equation solution for the 
rjNN system corresponding to different rj — N scattering amplitudes. Calculations have also been 
done using few body equations within a finite rank approximation (FRA) to generate T^d^rjd- 
The p — d final state interaction is included in the spirit of the Watson-Migdal prescription by 
multiplying the matrix element by the inverse of the Jost function. The rj — d interaction is found 
to be dominant in the region of small invariant rj — d mass, M^^. The p — d interaction enhances 
the cross section in the whole region of M^^, but is larger for large M^^. We find nearly isotropic 
angular distributions of the proton and the deuteron in the final state. All the above observations 
are in agreement with data. The production mechanism for the entire range of the existing data 
on the pd — > pdr] reaction seems to be dominated by the two-step model of rj production. 
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I. INTRODUCTION 



The current great interest in the ?7-nucleus interaction exists because of the attractive 
nature of the rj — N interaction in the s-wave and the consequent possibihty of the 
existence of quasi-bound, virtual or resonant ?7-nucleus states [2]. The exact nature of 
these states, of course, depends upon the precise knowledge of the rj — N scattering matrix 
at low energies. As the r/ is a highly unstable meson (lifetime ~ lO^^^s), this precise 
information is difficult to obtain directly. It can only be obtained from the eta producing 
reactions through the final state interaction. With this motivation, starting with the early 
experiments near threshold at Saclay on the pd ^ ^Her] and the pd pdrj reactions, 
measurements have been carried out near threshold and beyond at Jiilich and Uppsala using 
the COSY and Celsius rings respectively. In this series of experiments, the recent data 
on the pd pdrj reaction using the Wasa/Promice setup at the Celsius storage ring of 
the Svedberg laboratory, Uppsala are thematically complete and cover the excess energy, 

§, [Q = y/s — rriri — nip — rrid) ranging from around threshold to 107 MeV. The data 
] (integrated over other variables) include the invariant mass distribution over the whole 
excess energy range for the rj — d, rj — p and p — d systems and angular distributions for the 
proton, deuteron and the eta meson. Like the pd ^ ^He i] reaction, the (inclusive) rj — d 
invariant mass distribution exhibits a large enhancement near threshold and hence appears 
promising to study the rj — d interaction. The rj — p and p — d invariant mass distributions 
do not show any such enhancement. All observed angular distributions are nearly isotropic. 

Like in our earlier studies on the pd ^ ^He rj reaction our primary aim in this paper is 
to investigate the above mentioned data on the rj — d invariant mass distribution to obtain 
a better understanding of the rj — N a.s well as the rj — d interaction. We speculate, from 
our experience on the study of the pd ^ ^He rj reaction [4] , that in the region of low rj — d 
relative energy this set of data will be mainly determined by the t] — d interaction, though 
the three-body nature of the final state may introduce some uncertainty in this conclusion. 

We present a study of the pd pdrj reaction which includes the effect of the final state 
interaction. We have investigated two possible diagrams for the production mechanism: 
the direct mechanism and the two step process of t] production. The direct mechanism 
proceeds via an intermediate pn drj reaction with one of the nucleons in the deuteron 
as a spectator. The 1] meson in the two step model is produced in two steps, namely, pp ^ 
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diT^ and ti^ N 77 A^, hence involving the participation and sharing of the transferred 
momentum by three nucleons. The two step model for rj production was first used in and 
the data on the pd ^ ^He rj reaction was well explained. The vertices at the two steps have 
been described by the corresponding off-shell T-matrices. The T-matrix for 7r+ —>■ t] N is 
taken from a coupled channel calculation py], and that for pp dn^ is obtained from the 
SAID program provided by the authors of Ref. 6|. 

The final state interaction between the rj and the deuteron is explicitly incorporated 
through an 1] — d T-matrix, T^^. This T-matrix, which is required to be half-off-shell, is 
described in two ways. One choice involves taking a "factorized form" which is given by an 
off-shell form factor multiplied by an on-shell part given by an effective range expansion up to 
the fourth power in momentum. The parameters of this expansion have been taken from an 
existing recent relativistic Faddeev equation solution for the rjNN system \}\ corresponding 
to different rj — N scattering amplitudes. The off-shell form factor will be described in the 
next sections and is chosen to have a form without any adjustable parameters. The second 
prescription involves solving few body equations within the finite rank approximation (FRA) 
to obtain T^^. This approach has been used in literature for the f] — d, — ^He and — ^He 
systems [8| . We perform calculations for both the prescriptions using different models of the 
elementary coupled channel ?7-nucleon T-matrix which characterize them. 

The interaction between the t] meson and the proton in the final state, to a certain 
extent is contained implicitly in our calculations. This is due to the fact that we describe 
the 7r+ A^ — > 77 A^ vertex by a T-matrix, which has been modeled to include the t] — N 
interaction. This off-shell T-matrix treats the vr A^, riN and tt A channels in a coupled 

n 

channel formalism pj and reproduces the experimental data on this reaction very well. 

The effect oip — d final state interaction (FSI) is incorporated in the spirit of the Watson- 
Migdal FSI prescription jol , in which our model pd pdrj production amplitude is mul- 



tiplied by a factor which incorporates the FSI 



factor is taken to be the frequently used 
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3etween the proton and the deuteron. This 
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131] inverse Jost function, [J{p)] ^, where 



p is the relative p — d momentum. The assumption implicit in this approximation that the 
mechanism for the primary reaction be short ranged is very well fulfilled in the 77-production 
reactions. The momentum transfer in these reactions near threshold is around 700 MeV/c. 
We include FSI for both doublet (^5'i/2) and quadruplet (^5*3/2) p — d states. 

The ?7-nucleon T-matrix, which characterizes our calculations, is not precisely known. 
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Recent theoretical works on the np drj reaction [l^ conclude that the data on this 
reaction can be reproduced with the strength of the real part of the //-nucleon scattering 
ength ranging between 0.42 and 0.72 fm. In our earlier work on the pd ^ ^Het] reaction 
^, we found a good agreement with data, with the real part of the scattering length taken 
to be around 0.75 fm. This value was also found to be in agreement with the np ^ drj 
data in a i^-matrix calculation of the final state rj — d interaction in [iS^ . The same authors 
as in p^], recently performed a fit to a wide variety of data which includes the vrA^ — * vr A^, 
n N rjN, 'yN — > vrA^ and 'y N — > rj N reactions and gave their best fit value of 
the 77-nucleon scattering length, o^tv to be (0.91,0.27) fm [16|. The rj — d effective range 
parameters are given in [t^ for a^jv up to (1.07, 0.26) fm. Hence, in the present work we 
perform calculations with different models of the rj — N interaction, which correspond to 
three different values of the rj — N scattering length, ranging from arjN = (0.42 , 0.34) fm 
to (1.07, 0.26) fm. 

We find that the cross sections calculated using the two-step model and the above in- 
puts for the final state interaction reproduce most of the features of the experimental data 
reasonably well. ^ 

A theoretical effort to understand the Uppsala data [3] was made earlier by Tengblad 
et. al. |17j. In [l^] the contribution of three different diagrams, namely, the pick-up (a 
direct one-step mechanism of r] production), the impulse approximation and the two-step 
mechanism (here the 77 meson is produced in two steps via the pp d and N rj N 

reactions) to the cross section for the pd pdrj reaction is determined. The authors in 
\v\ conclude that the impulse approximation is in general negligible as compared to the 
other two diagrams, the two-step mechanism is dominant in the near threshold region and 
the contribution of the pick-up diagram (referred to as the direct mechanism in the present 
work) increases with energy and matches the two-step contribution at an excess energy of 
Q = 95 MeV. The latter conclusions regarding the contributions of the two step and pick up 
diagrams are in contrast to the findings of the present work as well as to existing literature on 
similar kind of reactions. We note here that the authors in 17|] do not include the final state 
interaction in their calculations in any way. They treat the kinematics and the dependence 
of the pion propagator (appearing in the two step model) on the Fermi momenta in an 
approximate way. The T-matrices which enter as an input to the two step model are simply 
extracted from experimental cross sections and are hence not proper off-shell T-matrices. 
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As a result of the above approximations, the authors in [17[ do not reproduce the observed 
enhancement in the rj — d invariant mass distribution near threshold, and unlike the observed 

isotropic distributions, find anisotropy in their calculated angular distributions. ^ ^ 

The contribution from the direct mechanism (or the so-called pick-up diagram of [17]) to 
the total cross sections is found to be about four orders of magnitude smaller than the two- 
step contribution at threshold in the present work. The one-step contribution does increase 
with energy (as also found in [13]), however, even at the highest energy for which data is 
available {Tp = 1096 MeV) it remains two orders of magnitude smaller than that due to the 
two-step model. This is in contrast to the observations in [l7|, where the two processes give 
comparable contributions at high energies. The difference of orders of magnitude between 
the two processes can be understood as a result of the large momentum transfer, q, in 
the one-step process. This q, which is very large in the threshold region (~ 840 MeV/c) 
continues to be large even at high energies. For example, it is ~ 600 MeV/c even at the 
highest beam energy of 1096 MeV. This finding of ours is very similar to the previous studies 
of the reactions involving high momentum transfer. For example, as mentioned above too, 
in [sl, for the pd Her] reaction up to 2.5 GeV beam energy, the authors comment that 
the one-step cross sections underestimate the data by more than two orders of magnitude. 
In yet another calculation [IS] of the cross section for the pd — s>^He X reaction (where 
X = 7] , 7]' , uj , (p) the two-step model was found to describe the data on these reactions up 
to 3 GeV quite well. In 19|], in connection with the pd ^HaK'^ reaction, the authors 



claim that for a beam energy of 1 - 3 GeV, the one-step mechanism predicts 2 to 3 orders of 
magnitude smaller cross sections as compared to the two-step mechanism. The cross sections 
obtained from the one-step model, in Ref. 17| are, however, reported to be only one order 
of magnitude less than those due to the two-step model at threshold and comparable to the 
two-step ones at high energies. 

In the next section, we describe the details of the formalism. In the subsequent sections 
we present and discuss the results and finally the conclusions. 
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II. THE FORMALISM 



The differential cross section for the pd pdr] reaction, in the center of mass, can be 
written as, 

da = — — — ^ dQp' \kp>\ dM^d \knd\ dQr,d 7: {\T\ ) (1) 

where ^/s is the total energy in the center of mass and kp and kp> are the proton momenta in 
the initial and final states respectively. M^^d denotes the invariant mass of the rj — d system 
and k^d and fi^^ denote, in the rj — d center of mass, the 77 momentum and its solid angle, 
respectively. Qp/ represents the solid angle of the outgoing proton. Angular brackets around 
|Tp in Eq. ([T]) represent the sum over the final and initial spins. 

The T-matrix, which includes the interaction between the rj and the deuteron is given by 

T = ( tpnd{kr,d), kp'] rup', rud' \ Tpd^pdr, \ kp, kd (= - kp)] rup, rud ) (2) 

where the spin projections for the proton and the deuteron in the initial and final states have 
been labeled as mp, m^, mp/, and nid' respectively. Tpd^pdrj is the production operator. 

The wave function of the interacting t] — d in the final state has been represented as 
iprid{krid)- tcrms of the elastic rj — d scattering T-matrix, T^^, it is written as 

//-I I , f {kr,d\Tr,d\q) / ^^ 

The second term here represents the scattered wave. It has two parts originating from 
the principal-value and the delta-function part of the propagator in the intermediate state. 
Physically they represent the off-shell and the on-shell scattering between the 77 and the 
deuteron. The on-shell part can be shown to be roughly proportional to the T] — d momentum 
and hence dominant at higher energies. The relative contribution of these terms in our case 
would be determined after we substitute the above expression for ipndiknd) in Eq- (12]) • We 
then get 

T = ( kr,d, kp'] TUp', rud' \ Tpd^pdr, I kp, kd{= -kp)] rup, uid ) (4) 

, \^ f dq { krid] rUd' | T^d | g; ^2' ) , ^ , „ , , 

+ y (27r)3 E{k d)- E{q)+ie ' ^^''^p' ' ^pd-*P<iv I ^rf! ^d) 

It can be seen that the T^^ here appears as a half-off- shell T-matrix. 
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FIG. 1: The two step process production mechanism for the pd — > pdrj reaction. 
A. The production mechanism 

For evaluating the 77 production T-matrix, ( | Tpd^pdr, I ), we assume a two-step mechanism 
as shown in Fig. 1. In this model, the incident proton produces a pion in the first step 
on interacting with one of the nucleons of the target deuteron. In the second step this 
pion produces an rj meson on interacting with the other nucleon. Both these nucleons 
are off-shell and have a momentum distribution given by the deuteron bound state wave 
function. To write the production matrix, we resort to certain standard approximations 

gram appearing in Fig. 1). The 
amplitude for the pN nd process, which in principal is off-shell, is taken at an on-shell 
energy. Considering the high proton beam energy, off-shell effects are not expected to be 



significant. The tt N n N process is inc 

n 

The production matrix is written as 



uded via an off-shell T-matrix. 

3, 



m's ^ ^ ^ TT TT 

where, the squared four momentum of the intermediate pion, = E^ — k^, with the energy, 
Et^, calculated at zero fermi momentum and k^^ = k^^ + kpi — k^jl ^ P. The summation is 
over internal spin projections and the matrix element {pn\d) represents the deuteron wave 
function in momentum space, which has been written using the Paris parametrization 21|. 
The factor 3/2 is a result of summing the diagrams with an intermediate 7r° and vr"*". 

The integral over the pion momentum in above includes the contribution from the pole 
as well as the principal value term. For the pion propagator itself, as we see, we have taken 
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-1/2 k -f-q 
p ^ 

FIG. 2: The direct process production mechanism for the pd pdij reaction. 

the plane wave propagator. This thus excludes any effect in our results due to medium mod- 
ification of this propagator due to other nucleons. This aspect may be worth investigating 
in future. 

The T-matrix for the intermediate pp — > Tr"*" c? process has been taken from an energy 
dependent partial wave analysis of the 7r+ pp reaction from threshold to 500 MeV |6|. 
The various observables in |6| are given in terms of amplitudes which are parametrized to fit 
the existing database. We refer the reader to jol and the references therein for the relevant 
expressions^f the helicity and partial wave amplitudes and the notation followed by the 
authors in 

For the n'^ n — > rjp sub-process, different forms of T-matrices are available. We use 
the T-matrix from [l| which treats the n N , i] N, and vr A channels in a coupled channel 
formalism. This T-matrix consists of the meson - N* vertices and the N* propagator as given 
below: 



with, 



-1 



where {a = tt , rj) are the self energies from the vrA^ and rj N loops. The parameters 

of this model are, qn* = 0.616, (3 = 2.36 fm~^ and Mq = 1608.1 MeV. This T-matrix 
reproduces the data on the Ti^n —>■ rjp reaction very well. 



Although the contribution of the direct mechanism (Fig. 2) is 

nn 

to the large momentum transfer involved in the process) [5|, 1 181 . 



mown to be small (owing 



13], for completeness, we 



calculate its contribution to the total cross section. The T-matrix for this mechanism can 
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be written as 



6 



{\Tpd^pd^\^) = -{\Tpn-.dr,{y/s^)\^) X |0d(g)P 



(7) 



where (pd represents the deuteron wave function in the initial state. The spin summed 
(I Tpn^drj P) is given in terms of the total cross section for the pn drj reaction by 



(TT{pn drf) 



2 rrip rrin nid 



TT S 



\Pi\ 4 



, (I Tpn^drjl ); 



(8) 



where Pi and p} are the initial and final momenta cm. system. The momentum transfer q, 
as shown in Fig. 2, is defined as, 

q = ^kp + kp> . (9) 



The total cross section, ax, for pn ^ drj reaction is taken from the experiments 22 1. 



B. Final state interaction 
1. 7] — d interaction 

This is incorporated through a half-off-shell r] — d T-matrix. We construct this T-matrix 
using the following two prescriptions: 

1. Factorized form of Tfjd 

In one ansatz we obtain it by multiplying the on-shell rj — d T-matrix by an off-shell 
extrapolation factor g{k' , k). Requiring that this T-matrix goes to its on-shell value 
in the case of on-shell momenta, we write 

T^.dik , E{ko) , k') = g{k , ko) T,.diE{ko)) g{k' , fco), (10) 

with g{p , g) — > 1 as p — > g. For a half-off-shell case, this obviously is the ratio of the 
half-off-shell to the on-shell scattering amplitude. 

For the on-shell rj — d T-matrix we use the effective range expansion of the scattering 
amplitude up to the fourth power in momentum, 

(11) 



Fik) = 4 + -Rk^ + Sk^ - ik 
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where F is related to T by 

T,d{k, k') = - I F,,{k, E{k), k'). (12) 

The effective range expansion parameters (A,R,S) are taken from a recent relativistic 
Faddeev equation (RFE) calculation of [3] ■ This calculation uses the relativistic version 
of the Faddeev equations for a three particle mNN system, where m is a meson and it 
can be an r], vr or a a meson. These particles interact pairwise, and these interactions 
are represented with separable potentials. The parameters of the rjN — ttN — aN 
potentials are fitted to the Su resonant amplitude and the n~p rjn cross sections. 
The rj — d effective range parameters obtained from these calculations are listed in 
for different sets of the meson-nucleon potentials. Each of these sets gives a specific 
value of the t] — N scattering length, which is also listed in [t]. 

Since the half-off-shell extrapolation factor g{k' , k^) is not known with any certainty, 
we choose the following two forms for it. 

(i) Following the method in [l5] for the final state interaction in the r] — d system, we 
express the off-shell form factor in terms of the deuteron form factor, 

g{k', ko) = j drj,{rk' l2)4>l{r)j,{rk,/2) (13) 

where for the deuteron wave function, ^^(r), we take the Paris parametrization. 

(ii) As a second choice, the form factor is taken to be the ratio of the off-shell rj — d 
T-matrix to its on-shell value, where both of them are calculated using the three 
body equations within FRA. The input to these calculations is the elementary t] — N 
scattering matrix, the details of which are given in the next Section. 

2. Few body equations within the finite rank approximation 

The other prescription oirj — d FSI involves the use of the half-off-shell rj — d T-matrix 
obtained by solving few body equations within the finite rank approximation (FRA). 
For the details of this formalism and the expression for the ?7-nucleus T-matrix, we refer 
the reader to our earlier works ^. To mention briefly, the FRA involves restricting 
the spectral decomposition of the nuclear Hamiltonian in the intermediate state to 
the ground state, neglecting thereby all excited and break-up channels of the nucleus. 
This is justified in the r] — ^He and possibly in the rj — ■^He case, but in //—deuteron 
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collisions, where the break-up energy is just 2.225 MeV, the applicability of the FRA 
may be limited. However, it should be noted that a comparative study [23] of the r] — d 
scattering lengths calculated using the FRA and the exact Alt-Grassberger-Sandhas 
(AGS) 2J] equations (which include these intermediate excitations) shows that they 
are not very different if the real part of the t] — N scattering length is restricted up to 
about 0.5 fm. 



2. p — d interaction 



We incorporate the p — d FSI in our calculations by multiplying our model T-matrix by 
the inverse Jost function, [J{p)]^^. We include the FSI in both the 1/2 and 3/2 spin states 
of p — d and restrict it to the s-wave. Since the p and d are charged we also include the 
Coulomb effects. Following standard procedure, we write the Jost function in terms of phase 
shifts and use the effective range expansion for the later. 

The complete expression for the s-wave inverse Jost function squared is written as. 



[Jo (V)] = [Jo {kpd)]Q + 



[(1 + ^)^o(M] 



-2 
D ■ 



(14) 



Here, to include the effect of the existence of one bound state, namely, the spin 1/2 state 
{^He), the doublet Jost function is multiplied by a factor (1 + ■^^^), where \Eb\ is the 
separation energy of ^He into p — d. Its value is taken to be 5.48 MeV. 

The expressions for spin quadruplet (Q) and doublet (D) [Jo (/cpd)]^^ are given by 



[Jo {kpd)]Q^ 

[Jo {kpd)]D^ 



4 



X 



X 



where. 



a = \ — 



3 Co fcprf 



q (^2 1.2 



ap 



1 + cot^ 6q 
1 

1 + C0t^5D 



(15) 
(16) 

(17) 



and af, and 6f, are defined as 

H' H' 



1 



2 7 kpd 



(18) 
(19) 
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where /i stands for either Q or D. The factor in above has its origin in the Con 
interaction. The phase shifts 5q^d are obtained from an effective-range expansion 25 . 



omb 



26|, 



C'^kpdCotS^ = + ^b^kl^ - (20) 



a,, 2 



7 = -TT— (21) 



hk. 



pd 

oo 2 

= E ( 2^ 2^ - ^^(7) - 0.57722 (23) 

n = 1 ^ ' 

Here rrired is the reduced mass in the p — d system, 7 the Coulomb parameter and a is the 
usual electromagnetic coupling constant. The values of the expansion coefficients a^, 6^ in 
Eq. f l20|) are taken as ag = 11.88 fm, Bq = 2.63 fm, an = 2.73 fm, and bo = 2.27 fm. They 
have been determined from a fit to the p — d elastic-scattering phase shifts in the relative 
p — d momentum range up to around 200 MeV/c 27 1. 

The above expression for the Jost function has the required property that for large p, 

Mp) ^ 1. 

iii. results and discussion 

Before we discuss the results of the present work, in order to highlight the FSI effects 
in the experimental rj — d invariant mass distribution we remove the phase space from the 
experimental da/dM^^d and plot in Fig. 3 the which is then given by, 

\f? = ■ ^ , (24) 

dMj^ci phase space 

where, 

phase space = — / dVlp' \kpi\ \krjd\ dil^^d (25) 

12 (27t)^ s\k„\ J 



12 {27lfs\kp^ 

as a function of the excess energy, Q^d = M^^d — fn-q — ^d, where M^^ is the invariant mass 
of the 7] — d system. In this figure we also show the plane wave result (i.e. Tpd^pdri does 
not include any FSI). The cross section, da/dM^d in Eq. fl2^ . is evaluated for each M^^d 
by performing an integral over the p — d centre of mass momenta, kpd- The range of the 
allowed values of kpd at each M^^ is shown by the hashed region. One clearly sees a large 
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FIG. 3: The ratio of experimental differential cross sections [3| to the phase space (Eq. (125]) ) as a 
function of the excess energy, Q^di along with range of p — d relative momenta, kp^ (hashed region), 
contributing to |/p at each Qr^d- 

enhancement in the experimental |/p near small values of Q^d^ which, most likely is due to 
the 7] — d FSI. We also observe a rise at large values of Q^d- Examining the range of p — d 
relative momenta which contribute to |/p at each Qrjd, one can see that this rise occurs at 
small values of kpd, indicating thereby the possibility of a large effect of p — d FSI in this 
region. 

In Fig. 4, we show two sets of the calculated |/p along with the experimental results 
for a beam energy of 1032 MeV. These results include only t] — d FSI. We limit the range 
of Qnd up to about 10 MeV, where, this effect is large. In Fig. 4(a) we show results for 
the factorized prescription with the off-shell factor generated from the deuteron form factor 
and the on-shell part arising from the relativistic Faddeev equation (RFE) calculation of 3] ■ 
The results are shown for three different sets of interaction parameters in the RFE. Since 
these sets give uniquely different values of the r] — N scattering lengths a^iN, we identify 
them by their corresponding o^at values. For the results presented here, these values are 
0.42 + iO.34 fm, 0.75 + iO.27 fm and 1.07 + iO.26 fm. We see that our results reproduce the 
enhancement seen in the experimental |/p at small values of Qr^d- The absolute magnitude 
depends upon the choice of the RFE parameters. It increases with arjN, which designate 
these parameter sets. The set corresponding to a^^r = 1-07 + iO.26 fm, gives results closest 
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FIG. 4: The calculated |/p along with the experimental results for a beam energy of 1032 MeV. 
(a) The results correspond to the factorized form of T^^ with the off-shell factor generated from the 
deuteron form factor, (b) The results correspond to T^d obtained from few body equations within 
the FRA. The data is the same as in Fig. 3. 



to the experimental values. 

In Fig. 4(b) we show |/p calculated using few body equations within the FRA, for t] — d 
7SI. These results are shown for three different inputs of the ri — N T- matrix taken from 
Iq . The choice of these T- matrices is such that their scattering length values are close to 
those used in Fig. 4(a). Though this model has the limitation of retaining the intermediate 
nucleus in its ground state in the r/— nucleus elastic scattering, the off-shell re-scattering 
effects have been properly included. If we compare Fig. 4(a) and 4(b), the two sets of results 
are similar. 

In order to check the sensitivity of the results to the off-shell form factor used in the 
factorized rj — d T- matrix, in Fig. 5(a), we show the |/p calculated using two different 
off-shell form factors. The on-shell Tj^d is obtained from RFE and the off-shell part is either 
treated with a deuteron form factor (solid line) or a few body FRA form factor (dash dotted 
line) as explained in section II B. The elementary r] — N T- matrix parameters required for 
the calculation of the FRA form factor are taken from the parametrization of Green and 
Wycech {l6|. Even though the results (as shown in Fig. 4(a)) corresponding to the a^^v = 
1.07 + iO.26 fm seem to be the closest to the data, to compare the effect of using different 
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FIG. 5: Comparison of the two form factors for OrjAf = 0.75 + iO.27 fm. (a) Effect of using two 
different off-shell extrapolation factors for rj — d FSI on (b) Two form factors as function of 

off-shell momentum (A;'). 

off-shell form factor, we choose the results corresponding to a^Tv = 0.75 + 10.27 fm. We 
make this choice such that we can compare the two calculations for the inputs corresponding 
to a similar rj — N scattering length. It should be expected then, that the off-shell form 
factors obtained from two different methods should not differ much. This is seen explicitly 
in Fig. 5(b) where the two form factors are shown as a function of off-shell momentum {k') 
for an on-shell value, /cq near the low energy peak in the rj — d invariant mass distribution 
(to be discussed in Fig. 9 later). 

Next, we include in our calculations the effect of the p — d FSI. This is done by multiplying 
the pd —>■ pdrj squared T- matrix (Eq. (jlj)) used above by the inverse Jost function squared 
in Eq. ( |T5i) and Eq. ( |T6i) . and integrating it over the allowed range (as shown in Fig. 3) 
of p — d momenta, kpd for each Qrid- We show these results in Fig. 6 for the RFE (with 
deuteron form factor) model oi r] — d FSI, for the parameter set corresponding to a^jv = 
1.07 + iO.26 fm. We find, that the p — d FSI affects the results in the whole region of Qr^d-, 
while the effect of — FSI is confined to small value of Qr^d- The large effect oi p — d FSI 
in the region of small Qr^d-, however, may not be taken with confidence as the value of kpd 
in this region is large (as shown in Fig. 3), where the s-wave effective range expansion (Eq. 
(!20|) ) for the calculation of Jost function might not be sufficient. In any case, it appears 
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FIG. 6: The proton-deuteron final state interaction effects on tlie pd pdr] reaction at the beam 
energy of 1032 MeV. The dashed hne shows the plane wave results and the dashed dot (solid) line 
shows the effect oi r] — d {rj — d p — d) FSI for Ujjn = 1-07 + iO.26 fm. 

that the effects of both the t] — d and the p — d FSI on the i] — d invariant mass distribution 
are significant. If we disregard the calculated p — d effect for small Qj^d, the rj — d and p — d 
FSI dominate in regions well separated from each other. 

Apart from the FSI, another important ingredient of our calculations is the two-step 
description of the production vertex. Because of the large momentum transfer, we believe, 
as has also been stressed in Ref. [iTI], that the angular distribution of the outgoing particles 
is probably more sensitive to the description of the production vertex. Inclusive angular 
distributions have been measured for all the three outgoing particles in the pd pdrj 
reaction. In Fig. 7, we show the calculated angular distributions for all the three outgoing 
particles along with the measured distributions. We show results without any FSI, with 
7] — d FSI and with t] — d and p — d FSI both included. As each angle has contribution from 
a range of Qr,d as well as kpd, the calculated results include integration of the cross section 
over these variables. We find that the observed nearly isotropic nature of the experimental 
angular distributions for the proton, deuteron and eta already gets reproduced by the plane 
wave calculations. The effect of both rj — d and p — d FSI is large and persists over all 
the angles. Their inclusion brings the magnitudes of the proton and deuteron angular 
distributions near to experiments. The magnitude of the eta distribution, however, does not 
seem to be affected much with the FSI. 

Experimental data also exist on the total cross section. In Fig. 8, we compare the total 
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FIG. 8: A comparison of the total cross section for the pd pdr] reaction calculated with the 
description of the production vertex as a two step mechanism and direct mechanism, along with 
the measured cross sections for a^^r = 1.07 + iO.26 fm 3, 28]. 
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FIG. 9: Contributions from the ofF-shell and the on-shell rj — d scattering in the final state. The 
results are for o^at = 1.07 + iO.26 fm with the inclusion of only the r] — d FSI. 

cross sections calculated including both the t] — d and p — d FSI with the measured cross 
sections. The results are shown with the factorized form of t] — d FSI with deuteron form 
factor for the set corresponding to rj — N scattering length equal to 1.07 + zO.26 fm. As we 
see, the calculated cross sections are in good agreement with the experimental data. 

In Fig. 8 we also give the cross sections calculated for the one-step direct mechanism 
(Fig. 2) mentioned in the previous section. Near threshold, these cross sections are about 
four orders of magnitude below those obtained from the two-step model and two orders 
of magnitude smaller in the high energy range. As mentioned in the Introduction, this 
observation is similar to that in other works involving large momentum transfer reactions 



19| . and is understandable because the momentum transfer continues to be large 



(~ 600 MeV/c) in the pd —>■ pdrj reaction even at an excess energy as large as 100 MeV. 

Now we make an observation about the importance of off-shell scattering in treating 
1] — d FSI near threshold. The scattering part of the rj — d wave function (Eq. ([3])), gets 
contributions from the off-shell as well as the on-shell scattering in the nucleus. To see 
quantitatively the relative importance of these two contributions to the cross section for the 
pd — > pdrj reaction, in Fig. 9 we show their contributions separately in the rj — d invariant 
mass distribution. These results include only the rj — d FSI generated from the factorized 
prescription using RFE and the deuteron form factor for the 1] — d T-matrix. We find that 
near threshold the off-shell scattering completely dominates the threshold enhancement. At 
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FIG. 10: A comparison of the calculated results including both the r] — d and p — d FSI with the 
experimental results. The results are for a^jv = 1-07 + iO.26 fm. 

higher excess energy, as expected, the on-shell contribution takes over. 

Finally we show the nature of agreement of our calculated results with the invariant r] — d 
mass distribution. In Fig. 10, we compare the calculated results including both the 1] — d 
and p — d FSI with the experimental results. The results are for a^iv = 1-07 + iO.26 fm 
calculated with the factorized prescription using RFE and the off-shell factor generated from 
the deuteron form factor. As we see the overall agreement is reasonably good. 



IV. SUMMARY 



The invariant rj — d mass distribution in the pd — > pdrj reaction has been studied 
by describing the production mechanism in terms of a two step model with a pion being 
produced in the intermediate state. The 7] — d final state interaction (FSI) has been included 
in (a) a factorized form involving an on-shell T^j^ and two types of off-shell form factors and 
(b) by solving few body equations within the FRA. The p — d FSI is included through a Jost 
function. The conclusions of this investigation can be summarized as: 

1. Experimentally observed large enhancement in the cross section near small rj — d excess 
energy, is reproduced by the rj — d FSI. The rise in the cross section at large Qr,d 
(which corresponds to a range of small momenta, kpd) can be accounted for by the 
p-d FSI. 
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2. Quantitative reproduction of the large enhancement requires t] — d FSI corresponding 
to large values of a^A?. In the present calculation it is around 1.07 + iO.26 fm. 

3. The calculations successfully reproduce the observed isotropic angular distribution 
of the proton and the deuteron in the final state. The total cross sections for the 
pd pdr] reaction are also well reproduced. 

4. The off-shell part of the rj — d scattering dominates near threshold. 

5. The results for two different choices of the off-shell extrapolation factor in the factorized 
form of the rj — d FSI are similar. 
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